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QCD thermodynamics II Frithjof Karsch 

1. Introduction 

The QCD vacuum has a complicated structure that manifests itself most prominently in the 
confinement of quarks and gluons and the spontaneous breaking of chiral symmetry. These non- 
perturbative properties of QCD are known to be temperature dependent and eventually will dis- 
appear at high temperature; at least in the limit of vanishing quark masses we expect that QCD 
undergoes a phase transition from a hadronic phase to a new phase of deconfined quarks and glu- 
ons in which chiral symmetry is restored. It has been speculated that there could be two separate 
transitions in QCD at which quarks and gluons deconfine and chiral symmetry gets restored ^ |2|], 
a scenario that has indeed been observed in lattice calculations of gauge theories with fermions 
in the adjoined representation [^]. In QCD with quarks being in the fundamental representation, 
however, it seems that at least for vanishing quark chemical potential there is a unique transition in 
the chiral limit at which quarks and gluons deconfine and chiral symmetry gets restored' . 

While in massless QCD the chiral condensate is a unique order parameter for chiral symmetry 
restoration, there is no counterpart for deconfinement. The Polyakov loop, which is an order pa- 
rameter for deconfinement in the limit of infinitely heavy quarks [^], is non-zero at all values of the 
temperature whenever quarks have finite masses. Nonetheless, the deconfining properties of the 
QCD transition are clearly reflected in the behavior of bulk thermodynamic observables, e.g. in the 
rapid rise of the energy or entropy density as well as in the sudden increase in fluctuations of light 
and strange quark numbers. The sudden change in the latter reflects the liberation of many light 
degrees of freedom, quarks and gluons, which dominate the properties of the thermal medium at 
high temperature. In the chiral limit these sudden changes go along with singularities in bulk ther- 
modynamic observables, the specific heat as well as quartic fluctuations of the light quark number 
diverge or develop a cusp. At the same temperature the chiral order parameter and its derivative 
with respect to the quark mass, the chiral susceptibility, show singular behavior. In this limit it 
is obvious that the singular behavior in observables related to deconfinement and chiral symmetry 
restoration, respectively, are closely related. To what extent this close relation persists also for 
non-zero values of the quark masses then becomes a quantitative question that should be answered 
through numerical calculations in lattice QCD. 

This became of particular interest in view of a recent calculation ||^ that suggested that there 
might be a large difference in the transition temperature related to deconfinement observables on 
the one hand and observables sensitive to chiral symmetry restoration on the other hand. It has 
been suggested that in the continuum limit this difference can be as large as 25 MeV. However, 
calculations performed with 0'{a^) improved staggered fermions [Q, ||] so far did not show such a 
large difference. 

In this write-up we will discuss some results from lattice calculations concerning the inter- 
play of deconfinement and chiral symmetry restoration. As these results have been presented in 
July/August of 2007 in a very similar format at the '4th International Workshop on Critical Point 
and Onset of Deconfinement' and at the 'XXV International Symposium on Lattice Field Theory' 
the write-up of these talks has been splitted into two parts. In the first part [^ we discussed recent 
results on the QCD equation of state. In this second part we will concentrate on results that can 



Large Nc arguments suggest that the situation could be more comphcated for non-zero quark chemical potential 
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give insight into properties of the QCD transition itself. We will focus here on a presentation of 
results obtained with &{a^) improved staggered fermion formulations. Results obtained with the 
1-link, stout smeared staggered fermion action [^ have been presented at both meetings separately 



[10] 



2. The transition temperature 

Before entering the discussion on deconfinement and chiral symmetry restoration, let us briefly 
summarize the current status of calculations of the QCD transition temperature using various dis- 
cretization schemes. The goal here is, of course, to determine the transition temperature in the 
continuum limit of lattice regularized QCD with its physical spectrum of two light and a heavier 
strange quark mass. While the heavier quarks, e.g. the charm quarks, may influence thermody- 
namics at high temperature [|n|, [l^, they are not expected to affect the transition temperature. In 
fact, even dynamical strange quark degrees of freedom seem to have little influence on the value 
of the transition temperature. Differences in the transition temperatures of 2, (2 + l)-flavor and 
3-flavor QCD still seem to be well within the current statistical and systematic uncertainty. This is 
in agreement with the observed weak dependence of the transition temperature on the light quark 
mass or, equivalently, on the light pseudo-scalar meson mass, 

roTcimps) - roT,{0) ~ A {rompsf , (2.1) 

with J ~ 1, A<0.05. To be specific we have used here the distance ro extracted from the static 
quark potential (see part I [^) to set the scale for T^. The weak quark mass dependence of Tc is 
consistently found in calculations with ff{a^) improved staggered fermions [0, ^^ 17] as well as 



with Wilson fermions [|13|, |14j, |15|]. This has been taken as an indication for the importance of a 
large number of rather heavy resonances for building up the critical conditions, e.g. a sufficiently 
large energy density, needed to deconfine the partonic degrees of freedom in QCD. As these heavy 
resonances are only weakly dependent on the quark mass values the light chiral sector of QCD 
may play a subdominant role for the quantitative value of the transition temperature, while it does, 
of course, control the universal properties of thermodynamic observables in the chiral limit. In 
2-flavor QCD for instance^, the scaling exponent d appearing in Eq. [2.1[ will be related to critical 
exponents (j3, 5) of 3 -dimensional, 0(4) symmetric spin models, d = 2/J55 = 1.08. 

In Fig. |](left) we show results on the quark mass dependence of the transition temperature 
obtained in calculations with the p4fat3 staggered fermion action on lattices with two different 
values of the cut-off, aT = 1/4 and 1/6 [^. As expected, in addition to the obvious quark mass 
dependence of Tc also a dependence on the cut-off, a, is clearly visible. Asymptotically the cut-off 
dependence is expected to be proportional to a^, i.e. we expect to find 

roTc{mps,N,) - roUO, -) ~ A {rompsY + B/N] . (2.2) 



^The discussion carries over to the light quark sector of (2+l)-flavor QCD. However, in this case one has to keep 
in mind that a second order transition point might be reached already at a non-zero value of the light quark mass. This 
effectively will change the universality class of the transition from 0(4) to Z(2) and modifies the singular structure of 
the free energy. 
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Figure 1: Left: Quark mass and cut-off dependence of the transition temperature calculated with the p4fat3 
staggered fermion action on lattices with temporal extent Nr — 4- and 6 [R]. Right: Transition temperatures 
determined in several recent studies of QCD thermodynamics. From top to bottom the first two data points 
show results obtained in simulations of 2-flavor QCD using clover improved Wilson fermions on lattices 
with temporal extent Nr = 8, 10 and 12 [ pj[ , Q and Nr — 4- and 6 Jl^, respectively. The remaining data 
points have been obtained in simulations of QCD with 2 light quark masses and a physical strange quark 
mass. They are based on calculations with staggered fermions using the asqtad action on Nr ~ 4-, 6 and 8 
lattices [Q, the p4fat3 action on Nt: = 4,6 [01 and 1-link, stout smeai-ed action on Nt = 4, 6, 8 and 10 lattices 
[ph. Circles indicate that the determination of the transition temperature is based on observables sensitive 
to chiral symmetry restoration, i.e. the chiral condensate and susceptibilities deduced from it. Squares 
indicate that observables sensitive to deconfinement have been used to determine the transition temperature, 
e.g. the Polyakov loop, its susceptibility and/or light and strange quark number susceptibilities. The diamond 
indicates that both sets of observables have been analyzed. With the exception of results presented in [ [l5| ] 
all calculations aimed at an extrapolation to the continuum limit (A^^ -^ °°) for physical values of the quark 
masses. All results have been rescaled to a common physical scale using ro = 0.469 fm | 



This ansatz generally is used to extrapolate to the continuum limit and to extract the transition 



temperature, Tc = Tc{0,oo). Of course, when using Eq. 12 for an extrapolation to the continuum 
hmit one has to make sure that the asymptotic scaling regime has been reached. In Ref. [^ the 
extrapolation is only based on two different values of the lattice cut-off, aT = 1/4 and 1/6, which 
may not be close enough to the continuum limit. This has been taken into account in the analysis 
performed in Ref. [^ by estimating a systematic error for the possible scaling violations. This lead 
to an estimate of the transition temperature T^. = 192(7) (4) MeV with the second error denoting an 
estimate for the systematic uncertainty in the extrapolation. An earlier analysis performed with the 
asqtad action on lattices with temporal extent Nt = 4,6 and 8 but smaller spatial volume, Na/N-c = 
2, lead to the estimate T^ = 169(12) (4) MeV [[I^. Both calculations currently get improved in a 
systematic comparison of simulations performed with the p4fat3 and asqtad action on lattices of 
size323 8[|8ll. 

In calculations with the p4fat3 action Polyakov loop and chiral susceptibilities have been ex- 
amined. The transition temperature has been determined by locating peaks in these susceptibilities 
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Figure 2: Disconnected part of the light quark chiral susceptibiHty and the Poly akov loop susceptibility (left) 
[^ and the quartic fluctuations of the light quark number (right) P3| ] calculated on lattices with temporal 
extent A^^ = 4 in simulations with the p4fat3 action. 



(see Fig.^left)). This lead to consistent results on larger volumes, although systematic deviations 
have been observed for smaller spatial volumes, e.g. for Na/N^ = 2. Calculations performed with 
the 1-link, stout smeared action [^ on lattices with temporal extent ranging from A^^ = 4 up to 
Nt = 10, on the other hand, suggest that there exist differences in the location of the peak posi- 
tions in chiral susceptibilities and inflection points determined for observables like the Polyakov 
loop and strange quark number susceptibility. These differences seem to become increasingly sig- 
nificant with increasing A^^. In this analysis a transition temperature related to chiral properties 
is determined to be Tc = 151 (3) (3) MeV while observables related to deconfinement suggest a 



transition temperature r^ ~ 175 MeV. For further discussion of these calculations see also [10]. 

In the Wilson formulation a discussion of chiral symmetry restoration becomes more involved 
than in the staggered case. For this reason only observables related to deconfinement, i.e. the 
Polyakov loop and its susceptibility, have generally been analyzed in studies with Wilson fermions. 
Compared to calculations of the QCD transition temperature performed with staggered fermions in 
(2 + l)-flavor [^ 0,[l^] and 3-flavor [ [l9| ] QCD, calculations with clover-improved Wilson fermions 
[13, 14, Ol] performed in 2-flavor QCD typically use rather large quark masses, mpsro>\. This 



makes an extrapolation to physical masses difficult. Nonetheless, a straightforward application of 
the scaling ansatz, Eq. |2.2[ for extrapolations of the transition temperature down to the physical 
values of the pion mass, mpsro ~ 0.32, yields results that are in reasonably good agreement with 
values determined within staggered fermion formulations. The current status of calculations of 
transition temperatures with Wilson and staggered fermions is summarized in Fig. [l](right). 



3. Deconfinement and chiral symmetry restoration 

The relation between deconfinement and chiral symmetry restoration in QCD has been dis- 
cussed since a long time [[Ij, ^. Although both phenomena seem to be related to physics on different 
length scales lattice calculations seem to suggest that both phenomena happen at approximately the 
same temperature even at finite, non-zero values of the quark masses when none of the symmetries 
related to confinement (Z(3) center symmetry) or chiral symmetry breaking {SUi{nf) x SUR{nf)) 
are realized exactly. In calculations with almost physical light quark masses and a physical value 



QCD thermodynamics II 



Frithjof Karsch 



of the strange quark mass this has recently been confirmed on lattices with temporal extent Nt = 4- 
and 6 in a detailed analysis performed with 0'{a^) improved staggered fermions [Q]. The gradual 
restoration of chiral symmetry with increasing temperature is signaled by changes in the light and 
strange quark chiral condensates and the corresponding susceptibilities, 



{VV)c, 



T dlnZ 



m„ 



Xn 



Td^lnZ 



m„ 



(3.1) 



where q = I, s for the hght and strange quark sector, respectively. The onset of deconfinement, on 
the other hand, can be examined through an analysis of the Polyakov loop L and its susceptibility, 

XL, 

^ ■lTrfl^(«c.H),o) > XL = Nl{{L')-{Lf), 



{ 



3M 



(3.2) 



"0 = 1 



where [/. ^^x g denote the gauge field variables defined on temporal links of a lattice of size N^N^. 
Some results for the light quark chiral susceptibility, Xmj, and the Polyakov loop susceptibility, Xl, 
calculated in a simulation with a physical value of the strange quark mass and a light quark mass 
that corresponds to a light pseudo-scalar mass of about 220 MeV {mi/nis = 0.1) on lattices of size 
16^4 are shown in Fig. ^left). 

As can be deduced from this figure the Polyakov loop susceptibility does not provide a strong 
signal for deconfinement in calculations with light dynamical quarks; the Polyakov loop itself is 
non-zero at all temperatures and rises smoothly through the transition region. This results only in 
a shallow peak in Xl, which nonetheless is in good agreement with the peak position in the light 
quark chiral susceptibility, Xm,i- 

In part I |^ we have presented results on bulk thermodynamic observables, e.g. the energy 
density (e/T^), as well as light and strange quark number susceptibilities (Xi.s/T^)- At least on 
lattices with temporal extent N^ = 4 and 6 they both rise rapidly in about the same temperature 
range. In both cases this is due to the deconfining nature of the QCD transition; as indicated in the 
introduction both quantities are sensitive to the liberation of many light quark and gluon degrees of 
freedom. In the chiral limit the sudden rise of, e.g. e/T^ and Xi.s/T'^ seems to be closely related 
to the singular behavior of the QCD partition function that arises from the restoration of chiral 
symmetry. To be specific let us discuss here the chiral limit of 2-flavor QCD (see footnote 2). The 
singular part of the free energy, fg, is controlled by a reduced 'temperature' t that is a function of 



temperature as well as the quark chemical potential Hq [[22||. The latter adds quadratically to the 
reduced temperature in order to respect charge symmetry at Hq = 0, 



MT,^q)=b-'Mtb 



l/(2-a)^ 



~f 
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with t 



T-T, 



+ c 



(3.3) 



where b is an arbitrary scale parameter and a denotes a critical exponent. As a consequence the 
specific heat as well as the quartic fluctuations of the light quark number, cf ~ {{^y) ~^{Ny))' 
show singular behavior at the critical point t = that is controlled by the same critical exponent a, 
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Figure 3: Preliminary results of the hotQCD collaboration Iq] for the strange quark number susceptibility 
calculated on lattices of size 32^8 using two different G{cp-) improved staggered fermion actions, asqtad and 
p4fat3 (left). The vertical lines indicate a band of temperatures, 185MeV <T < 195MeV, which charac- 
terizes the transition region in the A^^ = 8 calculations l|^] (see also Figs. 4 and 5). The right hand figure 
shows a comparison of calculations performed with the p4fat3 action on lattices of temporal extent Nx — 4, 
6 [|3|| and 8 |§1 and with the 1-link, stout smeared action for N^ = 4, 6, 8 and 10 [^. Note that different 
conventions have been used to define the temperature scale (see text) 



Unlike the Polyakov loop susceptibility the quartic fluctuations of the quark number thus provide 
a strong signal for deconfinement. This is shown in Fig. |2|(right). In fact, a comparison of C4 calcu- 
lated here in (2+ 1) -flavor QCD with light quark masses that correspond to a light pseudo-scalar 
(pion) mass of about 220 MeV [^] and earlier calculations in 2-flavor QCD with 10 times heavier 



quarks corresponding to a pion mass of about 770 MeV [ ]25| ] show that the quartic fluctuations rise 
strongly with decreasing quark mass. 

In Fig. 5 of part I [^ we have shown results for the strange quark number susceptibility, i.e. 
the fluctuations of strangeness number Xs ~ (^i^)' calculated with the p4fat3 and asqtad actions on 
lattices of temporal extent A^^ = 4, 6 ||T6|, ||] and 8 [§] in (2-1-1 )-flavor QCD and a light to strange 
quark mass ratio mi/rris = 0.1. The Nt = S results are preliminary results obtained by the hotQCD 
collaboration. They are discussed in more detail in [Q]. These calculations indicate a quite good 
agreement between results obtained with the two different 0'{a^) improved discretization schemes, 
although in particular at temperatures above the crossover region some differences show up. This 
is more clearly seen in Fig. ^left) where we compare the preliminary results obtained within both 
discretization schemes on A^^ = 8 lattices. These differences may be due to small differences in the 
choice of quark masses that define the constant line of physics along which the calculations have 
been performed and may partly be also due to differences in the discretization errors for both actions 
which may be about^ 6% for A^^: = 8. These differences as well as the cut-off dependence of results 
obtained on N^ = 4, 6 and 8 lattice with the asqtad and p4fat3 actions are, however, small when 



^In the infinite temperature limit deviations from the continuum value, Xffeem=0' ^^^ ^^ calculated analytically. 
For massless free staggered fermions on lattices with temporal extent A'x = 8 this yields Xfree,m=0 / xf^^e m=0 ~ ^'^^ 
(asqtad), 0.98 (p4fat3), 1.47 (l-link,stout). Like in the case of the pressure and other bulk thermodynamic observables, 
cut-off effects in the quark number susceptibility are 0[cr) improved for the asqtad and p4fat3 action and only start 
with 1 jN^ corrections in the infinite temperature limit. 
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Figure 4: The difference of light and strange quark chiral condensates normalized to its zero temperature 



value as defined in Eq. 3.5 (left) and the renormalized Polyakov loop expectation value (right). Shown are 
results from simulations on A^^ = 4 and 6 lattice obtained with the p4fat3 [El]] action as well as preliminary 
results for A^T = 8 obtained by the hotQCD Collaboration [|||. The upper axis shows the temperature in units 
of the distance ro extracted from the heavy quark potential. The lower temperature scale in units of MeV 
has been obtained from this using ro — 0.469 fm [pO|]. The vertical lines indicate a band of temperatures, 
185MeV <T< 195MeV, which characterizes the transition region in the Nx = S calculations. 



compared to results obtained with the 1-link, stout smeared staggered fermion action [Q] as shown 
in Fig. 0(right). The differences between the asqtad and p4fat3 calculations on the one hand and the 
1-link, stout smeared calculations on the other hand arise from two sources. For small values of A^^, 
the quark number susceptibilities calculated with 1-link staggered fermion actions overshoot the 
continuum Stefan-Boltzmann result at high temperatures and reflect the strong cut-off dependence 
of thermodynamic observables calculated with this action. This is well-known to happen in the 
infinite temperature, ideal gas limit and influences the behavior of thermodynamic observables in 
the high temperature phase of QCD (see footnote 3 and also Fig. 2 in [^). On the other hand 
the differences also arise from the different choice for the zero temperature observable used to set 
the temperature scale. While the temperature scale in the asqtad and p4fat3 calculations has been 
obtained from the static quark potential (the distance ro), the kaon decay constant has been used in 
calculations with the 1-link, stout smeared action. Of course, this should not make a difference after 
proper continuum extrapolations have been carried out. At finite values of the cut-off, however, one 
should make an effort to disentangle cut-off effects in thermodynamic observables from cut-off 
effects that only arise from a strong lattice spacing dependence in a zero temperature observable 
that is used to define a temperature scale. In this respect, the scale parameter ro extracted from the 
heavy quark potential is a safe quantity which is easy to determine; it has been studied in detail and 
its weak cut-off dependence is well controlled [^, 24]. 

Let us now turn our attention to observables sensitive to chiral symmetry restoration which. 



of course, is signaled by changes in the chiral condensate (Eq. |3.1| ). This also is reflected in pro- 
nounced peaks in the light quark chiral susceptibility as shown in Fig. ^ As the chiral condensate 
receives additive as well as multiplicative renormalization, one should look at appropriate combi- 
nations that eliminate the renormalization effects. An appropriate choice is to subtract a fraction of 
the strange quark condensate from the hght quark condensate and normalize the finite temperature 



8 
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difference with the corresponding zero temperature difference, 



. (w)/T-;|(wkr 



We note, that the strange quark contribution to this quantity will drop out in the chiral limit; A/ (7) = 
lim,„,^oA/,i(r) thus will become the standard order parameter for chiral symmetry restoration. 
This normalized difference of condensates obtained in calculations with the p4fat3 action on lattices 
with temporal extent N^ = 4, 6 [^ and 8 [§] is shown in Fig. 0(left). 



In the right hand part of this figure we show results for the renormalized Polyakov loop ob- 
tained in the same set of calculations. As can be seen the most rapid change in both quantities 
occurs in the same temperature range also on lattices with temporal extent A^^ = 8 [Ql- A cut-off 
dependence, which shifts the transition region to smaller temperatures, is visible in both observ- 
ables. It, however, seems to be small and correlated in both observables. 

The rapid change in the chiral condensate reflected in Fig. ^by the drop in A/ ^ (r), of course. 



is correlated to a peak in the light quark chiral susceptibility, Xm.h introduced in Eq. 3.1. This 
susceptibility actually is composed of two contributions, usually referred to as the connected and 
disconnected part, 

Ztot = Zm,l = Zdisc +Zcon (3-6) 

with 

Zoo. = -W{Df'{x,0)D^'{0,x)) . (3.7) 

Here Di denotes the staggered fermion matrix for the light quarks. In Fig. 5 we show results for 
the disconnected part of the light quark chiral susceptibility and the combined total susceptibility 
obtained in calculations performed on lattices of size 32^ 8 with the asqtad and p4fat3 actions. This 
is compared to the subtracted, normalized chiral condensate, Ai^s{T), calculated with both actions 
on the same size lattices. As can be seen the peak in Xh^.c/T'^ ^s well as Xiot/T^ obtained from 
calculations within both discretization schemes is in good agreement and corresponds well to the 
region of most rapid change in A/j(r). 



4. Conclusions 

Studies of the thermodynamics of (2-i-l)-flavor QCD with a physical value of the strange quark 
mass and almost physical values of the light quark masses have been performed at vanishing chem- 
ical potential with two versions of i?{a^) improved staggered fermions, the asqtad and p4fat3 
actions. Already on lattices with temporal extent Nt = 6 they yield a consistent description of bulk 
thermodynamics, e.g. of the temperature dependence of energy density and pressure. This also 
holds true for the structure of the transition region and is confirmed through calculations closer to 
the continuum limit performed on lattices with temporal extent N^ = i. These calculations yield 
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Figure 5: Subtracted finite temperature chiral conden- 
sates normalized by the corresponding zero tempera- 
ture quantity evaluated at the same value of the cut-off 
(top), the disconnected part of the light quark chiral 
susceptibility (middle) and the total light quark chi- 
ral susceptibility (bottom). All figures show preUm- 
inary results of the hotQCD Collaboration obtained 
with two different &(cp-) improved staggered fermion 
actions on lattices of size 11? 8 [H]. 

preliminary results for the transition temperature 
that may differ by a few MeV, depending on the 
observable used to identify the transition at non- 
zero quark mass values. In particular, on the 
A^T = 8 lattice no large differences in the determi- 
nation of the transition temperature arises from 
observables related to deconfinement and chiral 
symmetry restoration respectively. The prelimi- 
nary results of the hotQCD collaboration indicate 
that the crossover region for both deconfinement 
and chiral symmetry restoration lie in the range 
r = (185-195) MeV for A^^ = 8 and m//m, = 0.1. 
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